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SUMMARY 

The analytic description of an orbital navigation program that will 
simulate the estimation of position and velocity of an Apollo spacecraft 
while in orbit about the moon or the earth is presented. The inertial 
positions of a series of landmarks are observed and used as the naviga- 
tional information, which is processed by linear filter theory techniques. 
The landmark positions are also estimated and their uncertainties in 
location are reduced. 

The program is applied in this report to the problem of estimating 
the position and velocity of an Apollo spacecraft while in orbit about the 
moon. Using realistic initial values for the uncertainties in landmark 
positions, it is shown for a typical case that the RMS position and veloc- 
ity estimation errors at Lunar Excursion Module separation are about three 
times as high as the unrealistic, overoptimistic case of assuming perfectly 
known landmarks. The landmark estimation uncertainties are decreased by 
about two-thirds of their initial uncertainties. 


INTRODUCTION 

The uncertainty in the locations of lunar landmarks will affect the 
accuracy of spacecraft navigation while in lunar orbit. The effect of 
these uncertainties can be reduced by including the landmark vector into 
the system state vector. The resulting nine -component state vector 
includes the three components of spacecraft position, three of spacecraft 
velocity, and three components of landmark position. Linear filter theory 
(see Appendix A) is then used to obtain the best estimate of spacecraft 
position and velocity along with the best estimate of the landmark position. 
As the estimate of the landmark position improves, it approaches the actual 
position with corresponding improvement in spacecraft position and velocity 
estimates. The amount of improvement is limited by the accuracy of the 
optical instrument used in measuring the orientation of the landmark with 
respect to an inertial system with origin at the spacecraft. 

Such a scheme for navigation in the vicinity of the moon has two very 
useful functions: 

(1) The navigation function itself, which is of course, the 
immediate problem. 

(2) The mapping function, which enables lunar landmarks to be 
determined to a higher degree of accuracy. This means, of course, that on a 
subsequent mission, a set of landmarks that are known to greater precision 
will exist. 

The purpose of this note is to present an analytic description of a 
spacecraft orbital navigation program and a determination of the current 
estimate of Apollo lunar orbit navigational accuracy. 
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There are several error sources that are involved in lunar orbital 
navigation which should be considered. Numerical values assumed for these 
errors are listed in Table I along with the appropriate reference. These 
may be divided into three groups: 

(1) The landmark location errors 

(2) Navigator acquisition errors 

( 3 ) The navigation system errors. 

The landmark location errors may also be broken down into convenient 
classifications. The first source of error is the uncertainty in the loca- 
tion of the landmarks relative to the moon-fixed system. The second source 
is due to the uncertainty in the inertial location of the moon (earth-moon 
distance, the orientation of the moon on its spin-axis, and the uncertainty 
in the spin rate of the moon. 

The navigation system errors are due to the inaccuracy of the optical 
instrument even under ideal conditions, and to the mis-orientation of the 
inertial platform. For purposes of this study both of these errors have 
been lumped into the standard deviation of the instrument. 

The navigator acquisition errors should be considered alone, although 
for simulation purposes these errors will be included only in the landmark 
location errors. This follows to a large extent from the fact that these 
errors will likely affect the results most severely by not acquiring the 
landmark properly. These error's then show up as errors in landmark loca- 
tion. 


There are, of course, other sources of navigator errors to which it 
is difficult to assign numbers. For example, eye fatique, lack of sleep, 
and general mis-orientation of the astronaut due to his new environment; 
all would have effects on the actual navigation. However, these error 
sources are not considered in the present investigation. 


OPERATIONAL CONSTRAINTS # 

The operational constrains that must be considered in the simulation 
are visibility and sighting frequency. The only visibility constraint 
that will be imposed is that the landmark must be sunlit. The closer the 
moon is to being full, the larger is the number of landmarks available 
for navigation. Obviously, this constraint reflects itself back to launch 
date and flight time. 
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The opacecraft is assumed to be in an 80-n.mi. circular orbit. At 
this altitude the distance from a spacecraft to a landmark on the horizon 
is about 23 degrees [see sketch 1(a)]. In practice, lunar landmarks will 
probably be difficult to identify and acquire until the spacecraft is much 
closer to the landmark than 23 degrees. It will be assumed that no obser- 
vation will be made unless the landmark is within a five-degree cone of 
the spacecraft, as shown in sketch 1(a). For a five-degree cone it is 
assumed that a maximum of about four sightings can be made if the sightings 
are approximately one minute apart. 


THE SYSTEM EQUATIONS 

Tl.e system configuration is the position and velocity vectors of the 
spacecraft (trajectory state) plus the position vector of a single land- 
mark while the spacecraft is actually in a position to observe the land- 
mark (as defined in the preceding section). When the landmark is not 
visible the system will refer to the trajectory state only. 

It is assumed that a statistical correlation exists between the 
errors in the trajectory state and errors in the landmark position. 

However, when a landmark is no longer being observed, it is assumed that 
this correlation vanishes because of subsequent measurements on other 
landmarks . It is also assumed that the estimated state is close enough 
to the actual state so that the two states may be considered to be linearly 
related. The deviations of the measured quantities from their actual values 
actual values are linearly related to the deviation of the estimated system 
state from its actual value (Appendix B). With these last two assumpti: :is 
the system may be estimated by means of the optimal linear filter theory, 
(reference 2). It is further assumed that the correlation matrix of the 
errors in a lanmark's position is used to initiate the computation when 
the same landmark is next observed. 

The equations of motion for the actual state are: 


x(t) = £(x,t), with the initial conditions x(t Q ) + €_(t Q ) (1) 

and for the estimated state 

0 

x(t) = f(x, t), with the initial conditions x(t Q ) = x(t Q ) (2) 

where jc(t Q ) is specified and §^(t Q ) is a vector of random position and 

velocity deviations. In reality, £ is impossible to determine. For 
purposes of simulation € may be computed from the eigen-values and 
eigenvectors of the initial covariance matrix E(t Q ) and a set of normally 

distributed random numbers as shown in reference 3- 
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The vector x is called the state vector and may be expressed as 


x T = (r, v, x^) (3a) 


when a landmark is being observed; or as 


x T = (r, v) (3b) 


when no observations are being taken. 

The covariance matrix of the errors in the estimate of the state may 
be propagated by the matrix differential equation. 


E(t) = F(t) E(t) + E(t) F T (t), E(t Q ) specified 


(4) 


The covariance matrix is defined by, 

(5) 

( 6 ) 


where 


■ a 


(e e T ) 


e = 


A 

x - x 


From this definition it is seen that E is a symmetric matrix, which for 
this problem is 9 X 9 when a landmark is visible, 6x6 when no landmark is 
visible. For convenience E may be broken down into augmented matrices, 


E = 


*7 *VL 


T 

E 


VL E L 


(7) 
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Ey is the covariance matrix of the trajectory and is considered at all times. 

E^ is the covariance matrix of the landmark being observed, and is the 

matrix which represents the correlation between the errors in the estimate 
of the trajectory and those of the landmark. 

The matrix P is composed of the partial derivatives of f with respect 
to the state vector x. 

When no measurements are being made the state vector is updated by 
integration of equation (2) and the covariance matrix by equation (4). 
Between measurements on the same landmark on the same orbit, equation (4) 
is integrated in itj entirety; between measurements of different landmarks 
only is integrated. 


The optimal linear filter (Appendix A) is employed at a measurement time. 
The state is updated by equation (A8) and the covariance matrix by equation 
(A9). The initial conditions for (2) and (4) are then reset by these updated 
values and the integration of (2) and (4) proceeds until the next measurement. 

Although the program to be described refers to using lunar landmarks 
as navigational information, the analysis is equally valid fen* determining 
the orbit of a spacecraft about the earth. The coordinate system in which 
the state variables are described is as follows: the position and velocity 

of the spacecraft (r and v) are in a selenocentric Cartesian coordinate 
system. The landmark vector x_ is in a selenographic (moon-fixed) polar 
coordinate system. When x_ is expressed in its selenographic Cartesian compo- 
nents (xl ) it may be operated on by a transformation matrix A, that will 
express ^it in the selenocentric system, yielding r_ = 



E y expressing r_ in the selenocentric system. 


k - 0 (8) 

so that V and f^ of the matrix F in equation (4) are zero. f V f 2> 

and f^ are just the components of v; and f^, f and fg are just the 

gravitational accelerations due to the primary attracting body (the moon 
in this case) plus the perturbations due to other bodies and their 
nonhomogeneities . 
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THE LINEARIZATION OF THE SYSTEM EQUATIONS 

Let the actual trajectory be considered a nominal about which the 
equations of motion will be linearized. Taking first order deviations 
of equation (2), 


$£ r r it)s £a> 


( 9 ) 


where 

% x «)* 

( 10 ) 

From Appendix B, it is seen that the measured angles 9^ and which 

are the right ascension and declination of the landmark (sketch 1(b)) with 
respect to an inertial system with origin at the spacecraft, are not 
linearly related to the state vector but to first order, their deviations 
from their actual values are. That is, 

le - H S'* (id 


Considering only these first-order deviations, it can be stated that 
they are related after a measurement at time t^ as the truly linear variables 
in equation A8. 

That is, 


r x'u,)= sxu,) + K(*,)[sg(0- 


( 12 ) 



( 13 ) 





since s§ is the deviation of the computed value ^ from the true 
value 0. The deviation of the measured angle 0 from the true value is due 
to the noise in the measurement only. 



G- 6 


* OL 


(iM 


So (12) becomes 




£'(<> e £ <■*,) + KU,)[e-§ + !L\ ( 15 ) 


The matrix K is computed as in Appendix A. 
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THE MEASUREMENT ERRORS 

The measurement error vector £ may he computed by considering the 
magnitude of a normally distributed error vector J perpendicular to 
(appendix B) and also that this error vector is uniformly distributed 
between 0 and Ilf about ^ . This deviation may be expressed. 

d m f * 

^oc, d G$$ /f 
J*U 9 » d Si/t ^ 

where is normall y distributed with zero mean and y is uniformly distributed. 
The errors in the two angles are, 

- y- r m/f - (l6) 


OCj * * sin /? (IT) 


The elements of R (Appendix A) may now be determined from the 
measurement errors, , and , 




The distribution functions of V and are, 


/ 

y on - yzT<r t 


e - 


-o»< Y< — 


(18) 



M * aF , 


119) 


o & S* zv> 


Therefore, 

-2Zl 

^ dAf <0, 

also, 

«o3* - 57«, -<9 

which says that the errors are independent. 

ORBITAL NAVIGATION SIMULATION PROGRAM (NAV) 

Observing the form of equation (15) it is seen that the integrated 
values of the estimated state vector msy be used in the state estimation 
procedure rather than the deviations from a reference or nominal trajectory. 
Equations (2) and (4) are integrated up to the measurement time. At the 
measurement time the state vector and the associated covariance matrix are 
upciated by equations (15) and (A9). 

The calculated angles ft,, and ^ are computed with (Bl^) and (B15) from 
the estimated trajectory and the current best estimate of the landmark position. 
Equation (l) is integrated simultaneously with equation (2) to provide the 
actual trajectory, and assumed values of the actual landmark positions are 
input. These quantities are necessary in order to compute the actual angles 
© 1 and for use in the estimation equation,, equation (15) • 

In the orbital navigation program the covariance matrix E is replaced 
by a matrix W, where 



r . ww T 


(so) 
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This substitution is made to insure numerical stability. The covariance 
matrix may not remain positive definite after a large number of computa- 
tions. The matrix W is guaranteed to remain at least positive semi- 
definite. The entire estimation problem may be formulated in terms of the 
U -matrix as given in Appendix C, and reference 2. 


The covariance matrix is then computed in the program for information 
purposes from the W-matrix. The values of RMS position, velocity, and 
landmark position that may be given in this paper are computed from the 
trace of the covariance matrix. 


RMS position = 

RMS velocity = 

RMS landmark 
position 



( 21 ) 


The gravity model used in the simulation consists of the moon as the 
primary, or reference body; the perturbations due to moon's triaxiality; 
the earth; the second, third, and fourth harmonics in the earth's potential 
function; and the sun. The trajectories ere computed by the Encke integra- 
tion technique. 


SIMULATION RESULTS 

A computer simulation was made to estimate the position and velocity 
of the Apollo spacecraft, at approximately 80 n.mi. altitude circular orbit 
above the moon, from sightings on five lunar landmarks on the front side with 
3 sightings per landmark. Two rains were made for comparison. In the first 
rrn, figure 1, it »ras assumed that the landmarks were perfectly. In the 
second rain, figure 2, it was assumed that the landmarks were known only to 
the accuracies given in Table 1. Both runs use a .003 radian uncertainty 
for the combination of sextant and navigation system errors. 
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At four hears from Lunar Orbit Insertion, which is approximately the 
time for beginning of Lunar Excursion Module descent, the RMS position and 
velocity for the case of perfect landmarks (figure l) were 0.34 n.mi. and 
1.4 fps. At the same time, for the case where the landmarks were poorly 
Known, (figure 2) the RMS position and velocity were 0.9 n.mi. and 4.0 fps. 
Comparison of these results indicates the amount of error introduced if the 
landmarks are assumed to be known perfectly. 

It is seen from figure 2, that the RMS position and velocity increase 
between measurements and decrease when a measurement is made, as would be 
expected. The net effect of the measurements is a reduction in RMS position 
and velocity. From figure 2c it is seen that the IMS landmark uncertainty 
also decreases, and the reduction is from just below 8000 feet to less than 
3000 feet for the last three landmarks sighted. 

CONCLUDING REMARKS 

The analytic description of an orbital navigation program has been 
presented. This program will estimate the position and velocity of a 
spacecraft while in orbit about some central body using as navigational 
information the observed inertial positions of landmarks located on the 
central body. At the same time the program estimates the position of the 
landmarks with respect to a coordinate system fixed in the body. 

This program has been applied to the problem of determining the orbit 
of the Command and Service Module between the times of Lunar Orbit Insertion 
and Lunar Excursion Module separation . It is shown that if the landmarks 
are assumed to be perfectly known, overopuimistic estimates of RMS position 
and velocity of the CSM at LEI! separation will occur. Using realistic 
values of landmark uncertainties, the RMS position and velocity at LEM 
separation are about three tines as high as for the same case using 
perfectly known landmarks . The landmark RMS positions sure decreased by 
about a factor of two-thirds from their initial value. 
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APPENDIX A 


THE OPTIMAL LINEAR FILTER 
Given the linear dynamic system 

X - F(*)x. , t ^ < u > 

Where £(*<**> 5 C i.))= jU(o) is known and the 

£U £(«) ) * e it*) is given 

the system (Al) is observed through 

& = M(J )£ +■ « (A2) 

made at discrete times. Ot is white gaussian noise, with zero mean, 

£(sf ) -o 

also £Cglt)^ T (f)) - Slt-?)XU) 

and 

Kit) - 7e r Ct)>o 

The best estimate of the initial conditions for (Al) are ) 

which are given. The best estimate at a later time for x is. 


x a,) - £(■*•) 


(A3) 



A-2 


where, 

§u,t.) -ru) $(*,*»*, I 

The error in the estimate is defined as, 

e.tf) - xU)-jL(t) w 

The covariance matrix E is defined as, 

ect) • t(e.Ct)& T (.t)) £U)>o (A5) 

which is the solution to the Riccatti equation, 

EU)* E(t)EU) t ECt) E’it) < A6 > 

which is integrated with the initial conditions, 

EL*.) - £[s.ct.)e r u 0 )] ( at ) 

If a measurement is made at t^, the best estimate for is assumed 
to be, 

xU.)- XL*.)* /<(*,)[ §(*,)- &(*,)] (A8) 



A- 3 


where O ( 4 , ) is the measurement at t and 
measurement is expected to be. 




s what the 


The covariance matrix is updated at t. by 

r 1 l 


£ U.) ‘ £U.) -xa.) xu, )eu, ) 


(A9) 


The optimum value of the weighting matrix 


xu.) 


is. 


*U,)= £((.) X T U.)*~'U,) 


(A10) 


Uhere ’ MUi* XU, ) E(i,)H T a, ) +-XU) 


(All) 


K(t^) is optimum in the sense that it minimizes the trace of E(t^), 
i.e., it minimizes the sum of the mean squared errors in the estimate of x. 
After the measurement, equations (A3) and (a 6) are re-initialized and then 
integrated to the next measurement time 



APPENDIX B 
P“ORM OF THE H MATRIX 

The H matrix for this problem is a two-by-nine augmented matrix, 
since there are two measurements and the state vector is a nine vector. 

Specifically, 


H = 



(Bl) 


The vectors h t and /jt x are associated with the angles ^ and 


respectively. 

From sketch 1(b) it is seen that 


e t -- $(/) t-'.z 

Taking the first variation, 

s&i =(*&£) 

\ €>S ' 

But and, ^ 

SO . £ ss /? £ X ^ — $ Ar ( B^) 


Substituting BU into B3, 



( B2) 


(B3) 



B-2 


It is desirable to obtain 6 in selenographic latitude, yU , 
longitude, X, and altitude, •£ , rather than in cartesian form. 

Xu* C*m + *) **** 

i'll «. C <Jr *+A ) ear jtL 4th A 

X' it . ( Km*#) *l»yu 


Taking the variation, 

j -Xw ^ 

S Xu I . I - tit X 

£ tL*J I CRn-^A)OOt/t 

or in vector-matrix notation. 


tit 


t-l] "1 I S/U 

■t L l/Kmd I i X 


S u = 3 SX L 


Substitute b 6 into B5> 


(b6) 



(BT) 


(B8) 

(B9) 

(BIO) 


B7 becomes 

SOi * £/ €Aj + c[&A£ +di fjfd. 

- [ kl sf 4/] st - £ 



B-3 


since 





(HI1) 


which is the same form as (A2) except for the noise. To first order, the 
expected deviation from the true value of the measurement vector is linearly 
related to the deviation in the state vector. 

f $ - y S' x (HL2) 

whert. 



(B13) 


The only vectors in H that are really necessary 8re 0/ and 

since and are zero and and are related by an orthogonal 
transformation to b t and . Prom P8 it is seen that we need only the 
partial derivatives of ^ and with respect to the components. These 
are readily obtainable from the geometry, since from sketch 1, 


oos e, - A 


(BlU) 

7 > 


oos & x * 2- 


(B15) 



Therefore, it may be 

stated that. 


4 T = ( A 

- ( T* ' 

-£. . *> 
•p* 

( m6) 

i T x - ( , 

Ajfs., 2—) 
TS* S' ' 

(BIT) 


The orthogonality of the measurements ^ and is immediately 
verified from Bl6 and BIT hy noting that. 


A, 


(B18) 



APPENDIX C 


CONVERSION OF COVARIANCE MATRIX TO W-MATRIX FORMULATION 

Because of numerical inaccuracies, that is, round-off and truncation 
errors, the covariance matrix E may not remain positive definite after many 
computational cycles. It is convenient to replace E by a matrix W such 
that. 


E = WW T (Cl) 

The W matrix is guaranteed always to be at least positive semi -definite. 

The entire estimation problem may conveniently be formulated in terms of 
■che W matrix. 

The Riccatti equation (k) used to propagate the covariance matrix 
between measurements becomes. 


W = FW 

The covariance matrix is updated when a measurement is taken by 

T 

w* = ws 

where 

T -IT 

SS = I - ZM Z 1 

m m 

Z = W X H 
M = Z T Z + R 

The optimum weighting matrix K takes the form, 

K « WZM' 1 


(C2) 

(C3) 

(c4) 
(c 5) 

(C6) 

(Cl) 



TABLE I 


ERROR SOURCES FOR ONBOARD LUNAR 
ORBITAL NAVIGATION 

(All errors are 24 * values) 


1 . Landmark Location Errors: 

(a) Uncertainty i: moon-fixed ( selenographic ) system.* 

ARC Dist. Horizontal Vertical 

From ( 0 °, 0 °) Error Error 


30 ° 

60° 


259 meters 
670 meters 
1190 meters 


98O meters 
1000 meters 
820 meters 


(b) Landmark inertial orientation 

( 1 ) Moon axial orientation 

( 2 ) Moon rotation rate 

( 3 ) Earth Moon vector 


errors:** 

500 meters each 
in landmark latitude, 
longitude, and 
altitude. 


2 . Acquisition Error by the Navigator:*** 
1000 meters each in landmark 
latitude, longitude, altitude 


3 . Navigation System Errors: 

(a) Inaccuracy of sighting instrument:*** 

(b) Misalinement of IMU 

(c) Lack of target resolution 
Total error = 3 m rad (1^) 


*Reference 1 

**Reference 4 

***Assumed realistic values 
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Spacecraft orbit 



sketch la— Observational limits. 




RMS CSM position, n mi 
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3 sightings /landmark /orbit 

Figure la. - RMS CSM position uncertainty as a function of time from LOI - assuming no landmark errors. 
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Figure lb. - RMS CSM velocity uncertainty as a function of lime from LOI. - assuming no landmark error. 
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3 sightings /landmark /orbit 

Figure 2c. - Reduction in landmark uncertainties as a function of time from 101. 
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